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Abstract 

Effects of strong interactions on the two loop electroweak radiative corrections to the muon 
anomalous magnetic moment, = {g^ — 2)/2, are examined. Short-distance logs are shown to 
be unaffected. Computation of long-distance contributions is improved by use of an effective field 
theory approach that preserves the chiral properties of QCD and accounts for constraints from the 
operator product expansion. Small, previously neglected, two loop contributions, suppressed by 
a 1 — 4 sin^ 6w factor, are computed and the complete three loop leading short-distance logs are 
reevaluated. These refinements lead to a reduction in uncertainties and a slight shift in the total 
electroweak contribution to aj^^ = 154(1) (2) x 10^^^ where the first error corresponds to hadronic 
uncertainties and the second is primarily due to the allowed Higgs mass range. 
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I. INTRODUCTION 



Recently, experiment E821 at Brookhaven National Laboratory achieved an order of 
magnitude improvement (relative to the classic CERN experiments) in the determination of 
the muon anomalous magnetic moment, — [g^ — 2)/2 . The new world-average value for 
that fundamental quantity is [1] 

= 116 592 030(80) x 10"^^ . (1) 

Additional data currently being analyzed should further reduce the uncertainty. 

At the present level of precision, comparison with the theoretical prediction for from the 
Standard Model requires knowledge of hadronic vacuum polarization effects with accuracy 
of one percent. The most recent dispersion integral analysis [2] (see also [3]) based on data 
from electron-positron annihilation into hadrons and hadronic r decays demonstrates that 
the issue is unsettled: e+e" annihilation leads to a prediction lower by about 3 a than the 
value (1) while the prediction based on r data is lower only by 1.3 cr. Moreover, the vector 
spectral functions derived from e^e~ annihilation and from r decays differ significantly for 
energies beyond the p resonance peak (by more than 10% in some regions). It seems very 
difficult to explain such a large difference by isospin breaking effects. Thus, it appears that 
the data from e+e^ annihilation and from r decays are incompatible: so, no conclusion can 
be derived yet about a deviation from the Standard Model prediction. 

Since a real deviation from theory would signal the presence of "new physics", with 
supersymmetry the leading candidate, it is extremely important that all such hadronic un- 
certainties be thoroughly scrutinized and eliminated as much as possible before implications 
arc drawn. Toward that end, new e^e^ and r data from Frascati and the B factories will 
hopefully help to resolve this puzzling difference. 

Beyond the leading hadronic vacuum polarization effects, strong interaction uncertainties 
also enter via higher orders that involve quark loops. Quark loops appear in light-by-light 
scattering contributing in three loops as well as in two-loop electroweak corrections. The 
latter are the subject of this paper although the hadronic uncertainties there are certainly 
much smaller than those induced by light-by- light scattering. 

At the two loop electroweak level, hadronic uncertainties arise from two types of dia- 
grams, quark triangle diagrams related to the anomaly and hadronic photon- Z mixing. The 
first category has been previously studied in a free quark approximation and the more gen- 
eral operator product expansion. Although phenomenologically both approaches produce 
very close numbers they differ: particularly with regard to their explicit short distance de- 
pendence, i.e. \ogmz terms. Here, we show that this difference is due to an incomplete 
operator product analysis in the second approach. When corrected, unambiguous short- 
distance contributions result. We also take this opportunity to update the long distance and 
total electroweak contributions. 

In the case of photon-Z mixing, its two loop contribution to is suppressed by a factor 
l — Aiiiv? Oyv ~ 0.1; so, it is not as important. It can be evaluated either in the free quark 
approximation (sufficient for logarithmic accuracy) or via a dispersion relation using data 
from e^e~ annihilation into hadrons. The difference is shown to be numerically insignificant. 

Finally, having clarified the leading short-distance behavior of the two loop electroweak 
radiative corrections to a^, we can use the renormalization group to estimate higher order 
leading-log contributions which, due to an interesting cancellation, turn out to be very small. 

In the end, our analysis leads to a new, not very different, but more precise and better 
founded prediction for the electroweak contributions to a^. 
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II. ELECTROWEAK CONTRIBUTIONS TO a 



In the Standard Model (SM) the one-loop electroweak contributions to a^, illustrated in 
Fig. 1, were computed about 30 years ago [4-8]. They have the relatively simple form 



2 



(2) 



where Gn — 1.16637(1) x 10 ^ GeV ^ is the Fermi constant obtained from the muon lifetime 
and 9w is the weak mixing angle. 
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FIG. 1: One-loop electroweak contributions to 



For sin^ 9w we employ the on-shell renormalized definition 

2 

•2/1 _ 2 1 '^W /o\ 

where mz = 91.1875(21) GeV and the W mass is correlated with the Higgs scalar mass, 
niH, via loop corrections (for rrit = 174.3 GeV) [9] 

mw = ("80.373 - 0.05719 In 0.00898 In^ ) GeV (4) 

V 150 GeV 150 GeV; ^ ^ 

For rriH = 150 GeV, the central value employed in this paper, we must use mw — 80.373 
GcV (rather than the direct experimental value, rriyy = 80.451(33) GeV, which corresponds 
to a very small mn), for SM loop consistency. That implies 

= 0.2231 (5) 

and 

aj^(l-loop) = 194.8 X 10"^^ (6) 

The calculation of two loop electroweak contributions to was more recent and con- 
siderably more involved. It started with the observation by Kukhto et al. [10] that some 
two- loop electroweak diagrams were enhanced by large logs of the form ln[mz/mn). Those 
authors carried out detailed calculations for a number of such enhanced diagrams. They did 
not account, however, for closed quark loops. At about the same time, in Ref. [11] it was 
shown that for superheavy fermions, like the top quark, logarithms of their mass appear 
in corrections to magnetic moments due to triangle anomaly diagrams. Detailed studies of 
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all closed quark loops were included in the calculation of a|]^ in Refs. [12, 13]. Finally, in 
Ref. [14], the two- loop calculation of all logs as well as constant terms was completed. 

Two loop corrections to naturally divide into leading logs (LL), i.e. terms enhanced 
by a factor of \n{mz/mf) where m/ is a fermion mass scale much smaller than mz, and 
everything else, which we call non-leading logs (NLL). The 2-loop leading logs are (see 
Eq. (76) below) 



,EW/ 



a-" (2-loop) LL = 



24V27r2 

36 




1 + 



31 
215* 



4^2 ^2 



, mz 
m 



In 



mz 
mf 



(7) 



where a ~ 1/137.036, Nf = 3 for quarks and 1 for leptons, /| is the 3rd component of 
weak isospin and Q/ is electric charge. Electron and muon loops as well as non-fermionic 
loops produce the ln(mz/m^) terms in this expression (the first line) while the sum runs 
over F = T.u, d, s, c, b. The logarithm ln{mz/mf ) in the sum implies that the fermion mass 
m/ is larger than m^. For the light quarks, such as u and d, whose current masses are very 
small, rrif has a meaning of some effective hadronic mass scale. 

In Eq. (7) we have retained for completeness small contributions from the ^-Z mixing 
diagrams in Fig. 2 which were previously excluded because they are suppressed by (1 
for quarks and (1 — 4s^)^ for leptons. 
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FIG. 2: Contribution to from the 7-Z mixing induced by a fermion /. 



The more important terms in Eq. (7), those not suppressed by (1 — 4s^), were checked 
by Degrassi and Giudicc [15]. They used knowledge of well studied QCD corrections to 
6 — s> 57 decay and translated them into into QED corrections to via appropriate coupling 

changes. The only place that they erred was for the small — ^ (^Ij — 2QfS^,^ (1 — 4s^) In ^ 
terms coming from quark loops. Technically the difference is due to ^-Z mixing (Fig. 2) which 
is proportional to QfQ^ (electric charges of the loop fermion and the muon): in [15] it was 
given as Q/- 

In addition to leading logs, the NLL 2-loop contributions have also been computed [14]. 
They depend on known constants, the top quark mass (here taken to be 174.3 GeV), 
\ii{mz/mw) terms, and the as yet unknown Higgs mass, mn- For tjih — 150 GeV, those 
corrections are numerically given by [14] 



,EW, 



2-loop) 



NLL 



-6.0 ± 1.8 X 10~^\ 



(8) 



where the error allows mn to range from ttih — 114 GeV (its experimental lower bound) 
to about 250 GeV. We have included in Eq. (8) small NLL contributions, —0.2 x 10~^^, 
proportional to (1 — 4s^) mf/m^^ induced by the renormalization of the weak mixing angle. 
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FIG. 3: Effective Z^^* coupling induced by a fermion triangle, contributing to a 



When evaluating Eq. (7), one is confronted by the presence of hght u, d. and s quark 
masses in the logarithms. They were used to crudely regulate long distance loop contri- 
butions in Figs. 2 and 3, where QCD effects were ignored [13]. For m„^d = 300 MeV, 
nis — 500 MeV, rric — 1.5 GeV, and rrib = 4.5 GeV, one finds 

aJW(2-loop)ii = -(36.7 ± 2) x 10"^^ , (9) 

where the error is meant to roughly reflect low-momentum hadronic loop uncertainties for 
the u, d, and s quarks in Fig. 3. Together, Eqs. (8) and (9) provide the two-loop total 
electroweak correction a|^^(2-loop) = — 42.7(2)(1.8) x 10^^^, which together with Eq. (6) 
leads to the generally quoted Standard Model prediction [16], 

af^ = 152(4) X 10-^\ (10) 

where the error of ±4 x 10~^^ is meant to reflect the total uncertainties coming from hadronic 
loop effects, the unknown Higgs mass and uncalculated higher order (three-loop) contribu- 
tions. 

Refinements in the above analysis are possible on two fronts: improvement in the low- 
momentum contribution of hadronic loops and an estimation of the leading three-loop elec- 
troweak contribution (which is part of the overall uncertainty) . The Higgs mass uncertainty 
will be overcome with its discovery. 

The easiest hadronic loop improvement can be made in the quark contributions to 7-Z 
mixing pictured in Fig. 2. Those effects, embodied in the last part of the bracketed terms 
in Eq. (7) can be obtained via a dispersion relation using cr(e+e~ — > hadrons) data. Such 
an analysis has been performed for various low energy processes. It effectively leads to the 
replacement [17] 

-| E iV, - 2 s^) In ^ ^ -6.88 ± 0.50 (11) 

q=u,d,s,c,b "^1 

which is somewhat larger than the value (—5.95) obtained in logarithmic approximation 
with constituent m, and s quark masses. (It suggests that smaller quark masses might be 
more appropriate.) Because those 7-Z mixing effects are suppressed by 1— 4s^, the shift in 
from this modification is tiny, —0.02 x 10~^^ , and can be safely neglected. 
Low momentum loop effects in the light quark triangle diagrams of Fig. 3 are more impor- 
tant. It is clear that the use of effective quark mass as an infrared cut off is in contradiction 
with the chiral properties of QCD in the case of light quarks. Indeed, in the chiral limit the 
infrared singularity in the quark triangle does not go away: it matches the Goldstone pole 
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in hadronic terms. This refers to the anomalous part of the triangle, i.e. to the longitudinal 
part of the axial current. 

The issue of how to properly treat light quark triangle diagrams, was addressed originally 
in a study by Peris, Perrotet and de Rafael [12] within a low-energy effective field theory 
approach for tt, rj, and rj' mesons. More recently, Knecht, Peris, Perrotet and de Rafael [18] 
have reexamined the issue using an operator product expansion (OPE) and Ward identities 
as guidance. We find their approach to the anomaly related longitudinal part of quark 
triangles to be a valid improvement over the naive constituent quark mass cutoff of Eq. (7). 

Unfortunately, their rather sophisticated OPE analysis failed to properly address the 
short-distance behavior of nonanomalous, i.e. transversal, part of the light quark triangles 
in Fig. 3. In particular, they do not reproduce the complete Inm^ terms in Eq. (7). That 
difference was attributed to QCD damping effects in [18] which were claimed to eliminate 
the nonanomalous In mz light quark contributions. However, in our opinion, it points to a 
shortcoming in their analysis. 

In section III, wc address in some detail what modifications to the study in [18] are 
required to restore the proper short- distance behavior. We then employ the effective field 
theory approach to improve the evaluation of light quark diagrams in Fig. 3, thus refining 
their contribution to . 

Having verified the short-distance behavior of Eq. (3) we are also in a position to evaluate 
higher order leading logs via the renormalization group. That analysis is carried out in 
section IV, where the leading log 3-loop contributions to are determined. Such a study 
was previously undertaken by Degrassi and Giudice [15]. Although we find small differences 
with their analysis, in the end wc also obtain a very small leading log 3-loop contribution to 
. In fact, the result is consistent with zero, to our level of accuracy, due to an interesting 
cancellation between anomalous dimensions and beta function effects. 

Finally, in section V, we give a refined determination of for which the errors are 
reduced and the central value is slightly shifted due to improvements in our analysis. 



III. HADRONIC EFFECTS IN QUARK TRIANGLES 

An interesting subset of the two-loop contributions to are those containing fermionic 
triangles of quarks and leptons, see Fig. 3. The internal triangles define the one- loop Z*77* 
amplitude where the Z and one photon are virtual while the other photon is real and soft. 
Those same triangles produce the well-known anomaly part of the Z boson axial-current. 
For cancellation of the anomaly, one needs to sum over all fermions in a given generation. 

In the case of a^^, the cancellation between quarks and leptons in Fig. 3 is not complete, 
because of their different masses and interactions. In this section we give a detailed analysis 
of the general structure of the Z*^Y amphtude, paying particular attention to the effect of 
strong interactions on the quark diagrams. Perturbative QCD corrections to short-distance 
logs are shown to vanish and have an overall negligible effect for heavy quark diagrams. In 
the case of light quarks, an operator product expansion and effective field theory approach 
are used to improve the evaluation of long-distance QCD effects. These refinements lead to 
an update of the fermionic triangle loop contributions to . 
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A. Structure of the Z*77* interaction 



We begin our general analysis by introducing some definitions. The interaction La- 
grangian for the electromagnetic and Z-boson fields, and Zy^ is 



-^int — e>l'^ 2\i 



9 



4 cos 9w 



zi'y, j^^EQffU, j'-E^ifh.isf, (12) 



where Qf and Ij are the electric charge and the third component of weak isospin and we 
retain only the axial part of the Z-boson current. (The weak vector current contribution in 
Fig. 3 vanishes by Furry 's theorem.) 

The Z*77* amphtude T^^^ is defined as 



T,y^iJ d'xe^^^ {0\T{j,{x)£{0)}\im 



(13) 



That is equivalent to 

T..„ = e T,,y , T,,, ^-j d^x d'y e'^^-^^ (0| T{j^{x) j,{y) j,^(0)}|0) , (14) 



where e"' is the photon polarization vector. We consider the limit of small photon momentum 
k. The expansion of T^^, in k starts with linear terms and we neglect quadratic and higher 
powers of A;. In this approximation there are two Lorentz structures for T^^, consistent with 
electromagnetic current conservation, 



le 



47r2 



(15) 



The first structure is transversal with respect to the axial current index the second is 
longitudinal. 

The contribution of Z*77* to the muon anomalous magnetic moment in the unitary 
gauge where the Z- propagator is i{—g,j,v + (ln(lv / t^'z) / {(f — Tn\) can be written in terms of 



Aa 



EW 



a 

TT 

X 



d^g 



(27r)4 q^ + 2qp 



3 V q^mf^ , 



mr, 



wl - 



■Wt] + 



rrir. 



nir. 



■ Wt 



(16) 



Here p is the four-momentum of the external muon. For logarithmic estimates, a much 
simpler expression is sufficient. 



Aa 



EW ^ ^ '^l 
TT 87r2V2 



POO 

y^dg^ 



Wl + 



rrir 



rriz + Q'' 



Wt 



(17) 



where = —q^. Moreover, the same expression with the lower limit of integration set to 
zero works with a power accuracy (in m^/m^) in the case of a heavy fermion in the loop. 
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The one-loop results for Z*77* can be taken from the classic papers by Adler [19] and 
Rosenberg [20]. In Ref. [10] they were considerably simplified in the limit of small external 
photon momentum. One then finds the following one-loop expressions for invariant functions 
wl,t, 

"1 da a{l — a) 



^l-loop 



(1 - a)Q^ + m) 



2 ' 



(18) 



where the factor Nf accounts for colors in the case of quarks. 

We also independently calculated T^^ using Schwinger operator methods for the fermionic 
loop. It can be presented as the polarization operator describing the mixing of two currents, 

and il but with the fermion propagators taken in the external field with the constant 
field strength. In the fixed point gauge x'^A^ = this propagator has the form [21] 



+ 



(p2 



m 



2X1 



(19) 



Then straightforward calculations lead to the above expressions (18) for the invariant func- 
tions Wt,l{q^) ■ 

The corresponding two-loop contributions to were calculated in Refs. [12, 13]. Ac- 
cording to Eq. (16) one needs to integrate using the Wt,l given in Eq. (18). Let us consider 
the part wt,l[/] which is due to the loop of a given fermion / with the mass m/ at the range 
of 3> mj . The asymptotic behavior of u't,l[/] at large is 



1— loop 



[/] 



2«;^-'°°P 



[/] = 4/|iV^gJ 



1 



2m} 

In 



o 



rrii 



(20) 



At large we can use the simpler form (17) for integration. It is clear then that for 
the individual fermion loop the integral / dQ'^WL is divergent. This reflects the fact that 
the theory is inconsistent unless the condition of anomaly cancellation between leptons and 
quarks is fulfilled. This condition has the form 



f 



(21) 



within every generation. It means that at ^ mj the leading terms cancel out after 
summing over fermions and Wl ~ (ln(5^)/Q^ implying convergence for aj^^ . 

Note the difference between the Wl and wt parts in the integral (17). The first one does 
not depend on mz at all while for the second we have a cut off factor l/(Q^-|-m|). Therefore, 
the Wt part is never divergent, instead individual fermion loops produce log(m^/m/) terms 
in when mf <^ mz- On the other hand, the one-loop relation Wt = Wl/2 means the 
anomaly cancellation condition (21) leads to cancellation of the leading terms in wt 

as well. It results in absence of log mz terms when mf <^ mz for all fermions in the given 
generation. 



B. Hadronic corrections for quark triangles 

How good is the one-loop approximation for wl and Wt? This question pertains to strong 
interaction effects for quark loops. As characterized in [12] this issue brings about a new 
class of hadronic contributions to the muon anomalous magnetic moment. 
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Let us first discuss perturbative corrections to the ^77* amplitude at Q 3> rriq due to 
gluon exchange in quark loops. The longitudinal function wl is protected against these cor- 
rections by a nonrenormalization theorem [22] for the anomaly. What about the transversal 
function wt? If the ag corrections were present for wt then after summing over the fermion 
generation, the leading term would become a.s{Q)/Q'^ ~ log Q/Aqcd)- According to 

Eq. (17) this leads to terms in which are parametrically enhanced by log(logm2/AQCD)- 

It turns out, however, that the ag corrections in wt are also absent at Q ^ rUq due 
to the new nonrenormalization theorem proved in Ref. [23]. The proof, stimulated by the 
present study, is based on preservation of the relation wt — wl/'^ beyond the one loop, i.e. 
in the presence of QCD interactions. This relation holds only for the specific kinematics we 
consider in which the external photon momentum is vanishing. 

Our above discussion means that corrections are absent for both wl and wt in the 
chiral limit rUq — 0. When quarks are heavy the ag corrections can appear but with a 
suppression factor myQ"^ at Q ^ rUq. In application to it implies that perturbative ag 
corrections are absent for light quarks. For heavy quarks, the logarithmic terms which are 
due to Q ^ rUq are not renormalized but non-logarithmic terms regulated by as{mq) could 
appear due to the range of momenta Q ~ rUq. 

Next comes the question of nonperturbative corrections. For the heavy quarks these cor- 
rections, given by some power of Ag^p/m^ , are small. As discussed above, the perturbative 
strong interaction corrections governed by as{mq) are also small for heavy quarks. In partic- 
ular, this argument is applicable to the third generation, r, b and t loops, so the free quark 
computational results obtained in Refs. [12, 13, 18] are very much under theoretical control. 

The first and the second generations contain light quarks u, d, s for which the momentum 
range of Q spans the hadronic scale Aqcd where nonperturbative efi^ects are O (100%) and 
give unsuppressed contributions to aj^^. This problem has been addressed in the literature 
and two approaches were suggested. In Ref. [13] effective quark masses for hght quarks in 
one-loop expressions were introduced as a simple way to account for strong interactions. This 
mass plays the role of an infrared cutoff in the integral over Q. A more realistic approach 
to the relevant hadronic dynamics was worked out in Refs. [12, 18]. Unfortunately, some 
conceptual mistakes in applying the OPE (we are going to comment on them in more detail in 
Sec. Ill E) led to incorrect results. This is immediately obvious in the ultraviolet sensitivity 
of the results in Ref. [18]: the dependence on Inm^ is not suppressed for the first and the 
second generations where all the masses are much less than ruz- 

For light quarks nonperturbative corrections to Z77* are given by powers of ^cu/Q^ 
while perturbative ones are absent as we discussed above. Thus, in the range of Q of order 
mz the one-loop perturbative approach applies and suppression of the dependence on mz 
due to the cancellation of the log mz terms for between leptons and quarks in the first 
two generations is guaranteed. 

The actual interplay of nonperturbative effects for light quark contributions to rep- 
resents an interesting picture very different for the longitudinal wl and transversal wt parts. 
For the first generation, in the chiral limit {mu,d — 0) nonperturbative effects are absent in 
Wl and the l/Q"^ one-loop behavior in hadronic terms matches the massless pion pole. This 
is the 't Hooft matching condition, [24] as was pointed out in [12, 18]. However, nonper- 
turbative effects are crucial for wti where they are responsible for a transformation of the 
singularity at small Q into p and ai meson poles. 

The situation is similar but somewhat more cumbersome for the s quark in the second 
generation due to the U(l) anomaly (the r]' meson should be included together with 77 



9 



meson). Also chiral breaking effects due to nis are more important. 

Below we present a detailed discussion of perturbative and nonperturbative effects for 
different generations. 



C. Third generation effect for a 



EW 



As we discussed above the one- loop expressions in (18) work very well for the third 

generation where both perturbative and nonperturbative corrections due to strong interac- 
tions are small. Substituting wl^t from Eq. (18) into (17) we get for the sum of r, b and t 
contributions to aj^^ the following result [12, 13, 18] 



9 m? 



„ , mf 2 mi ( , mr 5 \ , , mi , mi 8 
-In^--^ (ln^ + - +41n^ + 31n-|-- 
.2 3y 3 



mr. 



, (22) 



where we neglected small corrections of order ?Ti^/m^ ;,^^^, m'^ i,/m%, and m%/m^ . Numeri- 
cally, 

G„, m^ 



a 



Aa^^[r,6,t] = — 



^ ^ . 30.3 = -8.21 • 10~^^ 



(23) 



which is properly included in the results of Eqs. (7) and (8). 

Following the discussion in Sec. IIIB, we estimate a perturbative uncertainty by adding 
terms of order of as{mg)/7i to logmg. It gives for the uncertainty 



a G^ml /16 ^ as{mt) as{mb) \ 



(24) 



where Ct and are numbers of order ±1. Using as{mz) — 0.11, we come to an estimate 



±- 



a 



87r2x/^ 



• 0.3 ^ ±0.1 • 10 



-11 



(25) 



D. First and second generations: logarithmic estimates 

In the case of the first generation, i.e. u and d quark loops together with the electron 
loop, the characteristic hadronic scales are provided by the p meson mass, = 770 MeV, 
(for the vector current) and by the ai meson mass, ma^ = 1260 MeV, (for the transversal 
part of the axial current). Therefore, for Q below m^, only the electron loop contributes to 

On the other hand, the longitudinal part of the axial current is dominated by the tt meson 
whose mass is small. This dominance means that the one-loop expression (18) for WL[e, u, d] 
(but not for WT[e,u,d\) works all the way down up to Q ~ 171^,. Thus, the contribution of 
WL[e,u,d] in aj^^ is strongly suppressed. 

Considering In (m^/m^) as a large parameter we see that with logarithmic accuracy the 
first generation gives for 

Aa^K u,d]^^^ dg^ = ^ In 4 = -1-08 • 10"- , (26) 
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where we do not differentiate between rup and rua^. 

Tlie case of tlie second generation, n,c,s, is more involved. The cancellation of fermion 
loops takes place at > 4m^; so, we take mj/^ = 3097 MeV as an upper limit of integration, 
and m<^ — 1019 MeV plays the role of nip for the strange quark. In the interval between 
and mj/^ the muon and s quark loops should be included in the integration over Q. In the 
interval between = 547 MeV and only the muon loop contributes to wjn but we need 
to account for the pseudogoldstone nature of the r] meson. In contrast to the first generation 
where the longitudinal part of the axial current had the same quantum numbers as the tt 
meson, we need to re-express the axial current of the strange quark as a combination of the 
SU(3) singlet and octet, 

1 - 1 - 

= -sjul^s = --{u'Ju^qu + (i7i/75C? + S7i/75s) + -^{u'yu^5U + d'y^-f^d - 2s-f^-f5s) . (27) 

The singlet part associated with rj' does not contribute to wl below m^/ (which is of the 
same order as m^) but the octet part does, so u, d and s loops should be taken with octet 
weights. In the last range, between m^^ and only the muon loop should be counted. 
Thus, we arrive at 

Aaf^\p,s,c] = ~ g^2^ {-fml ^^ i^Lifi^ + s]) 

= ^ 41n^^ + -ln^ + 31n^ = -5.64 • 10"^^ 28 

To go beyond the logarithmic approximation, we must account for higher powers of 
for the light quarks. The Operator Product Expansion (OPE), which we next discuss, 
provides a systematic approach for computing them. It leads to refined estimates of hadronic 
effects. 



E. OPE considerations 

At large Euclidean one can use the OPE for the T-product of electromagnetic and 
axial currents, 

%^^if d^xe^^^rO-^(x)i,^(0)} = ^4,„^...„^(?)Or-S (29) 

where the local operators are constructed from the light fields. A normalization 

point /X, which the operators and coefficients implicitly depend on, separates short distances 
(accounted for in the coefficients c*) and large distances (in matrix elements of Oi). The 
field is light if its mass is less than fi. In the problem under consideration this includes the 
electromagnetic field of the soft photon which can enter local operators in the form of its 
gauge invariant field strength F^/j = d^Ap — d^Aa- It also includes gluonic fields as well as 
lepton and quark fields (with masses less than fi). 
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The amplitude T^,^ is given by the matrix element of f^i, between the photon and vacuum 

T,. = (0| T,. |7(^)) = E^,.a,...aM (0| ^r-"' hik)) ■ (30) 

i 

Since our approximation retains only terms linear in k, the matrix elements are linear in 

fiiv = k^ei, — kye^. This means that only Oi transforming under Lorentz rotations as (1,0) 
and (0, 1) contribute. In other words, the contributing operators should have a pair of 
antisymmetric indices, O"^ = — Of". The amplitude T^j, is a pscudotcnsor so it is convenient 
to choose Of^ to be a pseudotensor as well, which is always possible using a convolution 
with e^iyjs- Moreover, the C-parity of O"'^ should be — 1 . Retaining only the contributing 
operators we can write the OPE expansion as 

^^.. = E {ct(5')(-?'(^;.+?m5'^(^;.-?.?'^(^;,) +cl(g') q^q^'Ol^} • (31) 

i 

Parametrizing the matrix elements as 

(o|of |7(fc)) = -^«,r^ (32) 

where the constants Ki depend on the normalization point we get for the invariant func- 
tions Wt,L 

WT,L{Q') = T.<^A<i')'^^- (33) 



1. The leading d = 2 operator 

Operators are ordered by their canonical dimensions di which define the inverse power of 
q in the OPE coefficients. The leading operators have minimal dimensions. In our problem 
the leading operator dimension is d = 2, 

Of^^^F^^^^^e-^^%As, (34) 

where the operator F°'^ is the dual of the electromagnetic field strength, and the numerical 
factor is chosen in such a way that = 1. Its OPE coefficients appear at one loop. Note, 
that this operator corresponds to the unit operator in the product of the three currents in 
(14). 

In considering the third generation, the normalization point can be chosen well below 
all the masses rriT^b/, such that the coefficients c^ j;^ are given by a one-loop calculation, 

Ct,l ['^, b,t]^ wt,l [t, h,t\ (35) 

where wl^tIt, 6, t] are given by Eq. (18) with summation only over / = r, 6, t. Indeed, 
taking matrix element of T^i^ we are back to the one-loop expression for the amplitude 
. Perturbative corrections are governed by Q;s(mb), nonperturbative ones are of order 
~ (Aqco/'^^b)'' (due to the operator FGG of dimension 6, see the discussion below). 

Such a choice of the normalization point is not possible for the light quarks u, d, s: to 
apply the perturbative analysis to OPE coefficients we should choose ii ^ Aqcd , i-e. the 
normalization point /i is certainly much higher than the light quark masses. On the other 
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hand, for Q ^ Aqcd we can choose ji <^ Q . For this range Aqcd h Q we can 
use perturbation theory to calculate the OPE coefficients. In particular, for the leading 
operator (34) the OPE coefficients lif] chiral limit m/ = are given by the one- 

loop expressions for = 0] (see Eq. (20) at m/ = 0). An interesting point here refers 

to dependence on fj, for Liirif — O]. This dependence is absent not only at the simple level 
of logarithmic corrections but also at the level of power corrections in fi'^/Q'^. 

To demonstrate this let us consider the one-loop calculation of T^,^ in the background 
field method. Using the propagator (19) at m = and doing the spinorial trace we come to 
the following expression 

/ ^4(f+g)2 [iP + ^)^P'Pf^ + (^^ + <i)^p'fA > (36) 

where we have omitted an unimportant overall factor. The integral over Euclidean virtual 
momentum p is well defined both in infrared and in ultraviolet domains (the logarithmic 
divergence at large p drops out because of antisymmetry of F^^,). It means that the dominant 
contribution comes from the range of p ~ since there is no other scale. A simple integration 
results in 

^[q,.q'Fp, + q,q''Fp^]. (37) 

To this one has to add the loop with the massive Pauli-Vilars regulators which is simple to 
calculate using the propagator (19). This contribution adds a polynomial term —F^i, to the 
expression (37) what restores the transversality for the electromagnetic current j^j^ and leads 
to the results (20) for Wt^l at mj = 0. 

The calculation above proves that up to power accuracy the OPE coefficients 
Cy^^[m/ = 0] = wt,l['^/ = 0] . The portion of the integral (36) which comes from the range 
\p\ < ji constitutes a correction of order jJ^ /Q^. However, even this correction is absent if 
the symmetry features of the theory arc preserved. Indeed, the polynomial part ^F^i, which 
came from the regulator loop is due to arbitrary short distances so it does not depend on 
H at all. The conservation of the electromagnetic current fixes the coefficient between the 
polynomial regulator part and the dispersive part (37). Separating the part of the integral 
(36) with \p\ < we break the conservation. Moreover, due to chiral symmetry there are no 
other kinds of corrections, perturbative or nonperturbative, to the one-loop result for the 
longitudinal coefficient c£[m/ = 0] . That property is completely valid when we deal with 
the flavor nonsinglet axial current, like U'jiy'j^u — dji^j^d in the flrst generation, which has 
no gluonic U(l) anomaly. For the invariant function WLlruf — 0] that feature follows from 
the Adler-Bardeen theorem [22] and 't Hooft matching condition [24]. In terms of the OPE 
coefficients it translates in the equality of c^[mj = 0] and WL[Tnf = 0] because all other OPE 
coefficients vanish at m/ = for the longitudinal part. 

For the transversal coefficient c^[mf — 0] the situation is more subtle. As shown in [23] the 
symmetry of the dispersive part (37) under the n, u permutation leads to the pcrturbatively 
exact relation 2c^[mf = 0] = c;^['^/ = 0] . Nonperturbatively this relation is broken at the 
level of ^qcd/Q^ terms as we will see from the OPE analysis. 

We can also account for corrections to c^ i^ due to fermion masses which break the chiral 
symmetry. They can be read off Eq. (20) for wt,l ■ The corrections are of the second order 
in m/ but logarithmically sensitive to the normalization point /i which replaces m/ under 
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the log in Eq. (20) when we translate to c^l, 




(38) 



Summation over /, say for the first generation e, u, d, leads to the lepton-quark cancellation 
of the leading terms and we must consider operators of higher dimensions. 

2. Operators of higher dimension 

The next operators, by dimension, are those of d = 3 

Of = -t fa^P^.f = 1 e-P^'qf a,s • (39) 

Chirality arguments show that their OPE coefficients contain mass m/ as a factor, so 
by dimension, oc nif/Q^. To calculate these coefficients it is sufficient to consider tree 
diagrams of the Compton scattering type. 



cl^24^'-^. (40) 



Taking matrix elements between the soft photon and vacuum states we produce the following 
terms in the invariant functions WT,L{q^)'- 

A('='^WL = 2 A^'^='^WT = E ^/ Qf^f^f ■ (41) 

^ / 

If we neglect effects of strong interactions, it is simple to calculate in one loop with log- 
arithmic accuracy using, e.g., the propagator (19) in the external field and the normalization 
point fj, as the UV regulator, 

i^f^-Qf^frufln^. (42) 
rrif 

Substituting this Kf in Eq. (41) we observe the full match with the l/Q^ term in Eq. (38), 
together the d = 2 and d = 3 operators reproduce the one-loop wl,t in Eq- (20). 

Note that this match is for the terms of second order in mass. In QCD due to spontaneous 
breaking of chiral symmetry the matrix elements of quark operators (39) are not vanishing 
at rUq = 0, instead they are proportional to the quark condensate {qq)^ = — (240MeV)^. 
The operators (39) played an important role in the analysis by loffe and Smilga of nucleon 
magnetic moments with QCD sum rules [25] . They determined by a sum rule fit the quantity 

V = = - (43) 

^ ^7r^Qg{qq)o (350±50MeV)2 ^ ^ 

dubbed as the quark condensate magnetic susceptibility. 

Actually, the OPE analysis together with the pion dominance in the longitudinal part 
leads to a relation for magnetic susceptibihty similar to the Gell-Mann-Oakes-Renner 
(GMOR) relation for the pion mass [26]. This relation derived in Ref. [23] has the form 

{ijiu + md) Kg = -ml Nc Qq (44) 
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The GMOR relation F^ml = 



— {rriu + Tnd){qq)Q allows us to rewrite (44) as 



N = -4^^Q.(^-^)oX, ^--4j^-- (335 MeV)^ - ^^'^ 

Although this value of % is in a good agreement with the QCD sum rule fit (43) its magni- 
tude is about two times higher than results of other approaches, see [23] for references and 
discussion. 

Let us go further by operator dimensions. Nothing new appears for d = 4: all operators 
of dimension 4 are reducible to the d = 3 operators due to the following relation, 

f(Du,-f^ - D^-fi_,)-f5f = -mffa^^-f5f . (46) 

For d = 5 wc have operators ffF°'^ and f'y^fF"'^ (with factors rrif again) and at = 6 
there are many operators of the type {fcr'^^'j^f) iff), F'^^TiGnyG^^ and so on. These 
d = 5, 6 operators produce terms in T^i,. A particular example is the following four- 
fermion operator which appears due to diagrams in Fig. 4, 

of = g7"75 Qtqqi^l'tq-{a^^), (47) 

where the quark field q has color and flavor indices, the matrices of color generators f^, 
a = 1, ... ,8 act on color indices, and the weak isospin and electric charge Q are the 
diagonal matrices in the flavor space. This operator enters in the OPE (30) with the following 




coefficients 

^_16™^_ 4 = 0. (48) 

Note that our consideration of four-fermion operators is similar to Ref. [18]. Note also 
that the operator Oq contributes only to the transversal function wt — consistent with 
absence of nonperturbative corrections to in the chiral limit. Moreover, in the chiral 
limit, the d = 6 operators are next-to-leading after the leading d = 2 operator, showing that 
parametrically the leading nonperturbative corrections are of order Aq^j^/ . 

The matrix element of (47) between the vacuum and the photon states can be found 
assuming factorization in terms of the quark condensate (gg),, and the magnetic susceptibility 
Kg given in Eq. (45). It results in the following piece in wt for the light quarks in the first 
{u, d) and second (s) generations, 

Afrf-fi^ r ,1 oA(rf-fi^ r^ ?>2t: a s {qq)l (0.71 GeV)'' 
A('^-^)«;t[«, d] = -SA^'-'^wtIs] = ^ = -a, ^ . (49) 
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We can use this as an estimate for the terms in wt neglecting the FGG operators 

which enter with smaller coefficients (they appear in one loop while Oq is due to tree level 
diagrams). We also neglect the anomalous dimension of Oq. However, we have in mind 
that this anomalous dimension is rather large and positive and considerably compensates 
the running of ag which therefore can be taken close to 1 for estimates. 

Summarizing the consequences of OPE for the u, d and s quark loops in the chiral limit 
we get 

Wt[u, d]m„^j=o 

The longitudinal part given by the leading d — 2 operator F^^, has neither perturbative ag 
corrections nor nonpcrturbative ones, and the pole matches the massless pion. So the 
cancellation of w^' and wf^ in the ffist generation is exact in the chiral limit rriu^d = "^e = 0. 

As we discussed above the leading operator contribution to the transversal part has 
no perturbative corrections either [23]; however, nonpcrturbative corrections are present. 
Their signal in the chiral limit is very clean: the lowest masses in the vector and axial vector 
channels are nonvanishing in contrast with the pion in the longitudinal part. In Sec. Ill F 
below we present a resonance model for wt[u, d] consistent with the OPE constraints. Thus, 
there is no complete cancellation in the sum of w^' and although this sum decreases as 
at large Q. 



n r 1 1 (0.71 GeV)4 l\ 



3. Comparison with the OPE analysis in Ref. [18] 

It is convenient at this point to discuss a comparison of our approach with the analysis 
of the M, d quark loops of Ref. [18]. In essence, it is claimed there that the leading large Q 
behavior of the transversal part wt[u, d\ (and wt[s] as well) is 

■^t[u, 0^ ^ (51) 

in the chiral limit in contrast with the l/Q"^ perturbative behavior. As a result, the l/Q"^ part 
of wt[u, d] is absent while the l/Q"^ part of WT[e\ is present, so the quark-lepton cancellation 
is destroyed and, consequently, spurious logm^ terms appear in a}^ . 

To pinpoint the origin for such a dramatic difference between our approaches let us 
notice first that in [18] the authors consider T^^^,, the vacuum average of the product of 
three currents defined in Eq. (14), as a primary object. In this approach they do not have 
the electromagnetic field entering into local operators and our leading operator F^p does not 
appear. However, they have to consider the OPE for the product of three currents instead 
of two. For Tfj^^i, one can derive the following expansion: 

T,,. = (0 1 Ec;7.a,....,(?,^) 0r-"'(o) + ^ / ^'y^-'''T{f,M3,{y)] I 0) . (52) 

i •' 

The first part, which contains the local operators, accounts for emission of the soft photon 
from short distances. The second bilocal part, which contains the operator T^,^ defined in 
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Eq. (29) and the soft momentum current j^, accounts for the soft photon emission from large 
distances. Our relation (30) conveniently includes both parts in the local form: the first one 
is due to operators Oi containing the electromagnetic field strength Fa/s explicitly and the 
second is due to the operators without Fq,^ . Moreover, our leading operator Fq,^ corresponds 
to the unit operator in the first, local, part on the r.h.s. of Eq. (52). 

Once we have established this correspondence it is simple to see what is missing in the 
analysis of [18]: they did not account for the local part in Eq. (52). It is the unit operator in 
this part {F^p in our formalism) which gives the leading contribution and its coefficient 
follows from the perturbative triangle. This corresponds to soft photon emission from short 
distances as we discussed above. 

Note, that the unit operator in the local part of Eq. (52) is leading for both, longitudinal 
and transversal structures in the amplitude T^-yy, so its omission should result in an error 
for the longitudinal contribution as well. This did not happen in [18] because they did not 
apply the OPE analysis to the longitudinal part of the amplitude, instead fixing it by the 
anomaly. A comparison with the expansion (52) is particularly simple for the longitudinal 
part at mj = 0: the second bilocal term containing T^i^ vanishes and the only surviving 
operator is the unit operator in the first term. 

All the OPE subtleties discussed above should not screen a conceptually very simple 
situation: the short distance behavior given by free quark loops should not be changed in 
QCD by large distance effects. 



F. First generation 



We are now well prepared to calculate the contribution of the first generation to with 
an accurate account of hadronic effects. The invariant functions wt,l for the first generation 
is the sum of tyr,L[e], 'Ji'T,L[ii] and WT,L[d] ■ For the electron 



WL[e] =2wT[e\ = - — 



(53) 



where we neglected the electron mass. Expansions at large in the chiral limit for Wt^l[u, d] 
are given in Eq. (50). Hadronic effects modify wt,l[u, d] . Modifications are minimal for the 
longitudinal function wl[u, d]: the position of the pole is shifted to due to the explicit 
breaking of the chiral symmetry by quark masses,^ 

To find the contribution of WL[u,d\ to one needs to use the more accurate Eq. (16), 
rather than Eq. (17), because the integral is dominated by momenta Q ~ ttIt^ comparable 
with m^, 

i\a,Ae,u,d\ — 

"2 



/ da{l + a)\nA 
Jo 



+4 



m 



m 



1 1 ffi^ 

dail-aVlnA- - \n ^ 
3 



ml ^ 9 



(55) 



^ It is just this shift which allows one to derive [23] the expression in (44) by comparison of the terms 
with the OPE. 
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where ^4 = a + (1 — aY{m^^/rn^). Numerically it gives 

Aa^[e, u, d] = -- • 3.53 = -0.96 • 10"^^ . (56) 

The transversal function wt[u, d] can be modeled as a linear combination of two pole 
terms: one is due to the p(770) vector meson, another due to the ai(1260) axial vector 
meson, 

1 



m\ - ml ml- ml 



(57) 



"p ' ai. 

The residues in this expression are fixed by two conditions at large Q which follow from the 
OPE expression (50) plus the d — Z terms (41) breaking chiral symmetry. The first condition 
is on the coefficient of the leading 1/Q^ term, the second condition is for the coefficient of 
l/Q^. The term XjQ^ in (50) allows for an extra test of the model. The expression (57) 
gives — (0.96GeV)'' to be compared with — (0.71GeV)^ in the OPE based (50). Agreement 
is not extremely good but the right sign and order of magnitude are encouraging. Since the 
OPE \jQ^ estimate is very approximate, we use (57) for numerical estimates. 

For the integral over Q defining the contribution of u'ri'ii, o?], we can use the simpler 
expression (17) neglecting m^/m^ corrections, 

AaJle,»,d] = -^^|l.4-^^ln44|. (68) 
which gives numerically 

Aa^[e, u, d\ = -- ■ 4.88 = -1.32 • 10"^^ . (59) 



Overall, the first generation contributes to of^ 



EW 



Aa^W[e, d\ = -- ■ 8.41 = -2.28 • 10"^^ , (60) 

which is to be compared with the constituent quark model result [13], 



EWr Jl _ " G^,m\ 



^1 ^]free quarks 



n 87r2V2 



ml 17' 
In— ^ + — 



m^.m'i 



-4.0 • 10-^^ . (61) 



The refined result is about 1/2 of the constituent quark model value. It represents our main 
phenomenological finding. The primary reason for the shift is a deeper extension into the 
infrared due to quark-hadron duality for the longitudinal function WL[u,d]. It leads to a 
stronger quark- lepton cancellation for - the effect noted in Ref. [18]. 

What is the accuracy of the result (60)? Most of the model dependence is related to 
the description of the transversal function Wx[u,d]. For the longitudinal contribution, the 
analysis is rather solid. To get an idea of the accuracy, we consider variations when the p and 
Oi masses are in the intervals 500 — 1000 MeV and 900 — 2000 MeV. We found that deviations 
from the result (60) arc within 10%, i.e. of order of ±0.2 • 10~^^ . This high level of stability 
is related to the fact that the main contribution to Aa^[e, m, rf] in Eq. (58) comes from the 
unambiguous logarithmic term ln{m'j^/m'j^); it gives —1.08 ■ 10^^^ out of —1.32 ■ 10^^^ . 

So, in total, this analysis increases by 2 x 10~^^ relative to the free quark calculation 
and significantly improves its reliability. 
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G. Second generation 



The second generation contains both hght, s, and heavy, c, quarks which should be 
treated differently. For the light s quark we use the approach similar to the case of u, d 
quarks. For the longitudinal function wl[s] , as it was explained above in Sec. HID, one 
must include both singlet, rj' (960) , and octet, 7^(550), pseudoscalar mesons. 



For the transversal function the model is 



wt[s] = -- 



3 ml 



m". 



ml 



mz 



[Q^ + ml g2 + ml 



(62) 



(63) 



where 0(1019) and /i(1426) are isoscalar vector and axial vector mesons relevant to the ss 
channel. Integrating wl^t[s] and adding the known expression for the c quark and muon 
contribution we get for the second generation 



Aa^^[/.,.,c] = -- 



a Gf^mf^ 



TT Stt^v^ 13 



?ln^ 



?ln!^ 



^ 1 ml - ml ml 



mz 



m, 



3 ml 



mz, 



mz 



+ 41n^ + 31n^ - — 
ml mf, 9 




Numerically it constitutes (at mc = 1.5 GeV) 



Aa^«[/x,.,c] = - 



a G^mf^ 



17.1 = -4.63 • 10' 



-11 



TT 87r2v^ 

This numerical value practically coincides with the free quark calculation [13], 

47 Stt^ 



[/X, c, s] 



free quarks — 



a G^m^ 



In 



mz. 



m^m^ 



5+6 



= -4.65 • lo- 



ll 



(64) 



(65) 



(66) 



Two reasons for such a good agreement. First is the smallness of the strange quark contri- 
bution - its electric charge is smaller - so, hadronic details are not so important. Second 

is that the effect of cancellation between leptons and hadrons in the longitudinal invariant 
function w^is much less pronounced that in the first generation because of larger masses of 
T] and Tj' . 

The result is more sensitive to the c quark parameters. If we take 1.3 GeV for its mass 
instead of 1.5 we get (—4.32) • 10~^^, i.e. a change by 0.3 • lO"^-*^. Another source of the QCD 
corrections for the heavy quarks is perturbative gluon exchanges in the quark triangles. This 
estimation is similar to to one we did in Sec. IIIC, we substitute log mc by as{mc)/TT, 



a Gf,ml 
TT 87rV2 



8 



TT 



0.2 • 10-^^ , 



(67) 



where we used ^^(mc) ~ 0.3 . 

We conclude that the uncertainty coming from the second generation is small, about 
±0.3 X 10^^^, and related mainly to charm quark parameters. Overall, the total hadronic 
loop uncertainties in are well accounted for by an error of ±1 x 10~^^. 
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IV. LEADING LOGARITHMS: RENORMALIZATION GROUP ANALYSIS 



It was pointed out in [14] that once the leading log short-distance two-loop corrections 
to of order Inimw /tn^) are completely known, a renormalization group (RG) analysis 
can provide all leading log terms of the form [a\n{m\^ / ■> ^ = 2,3, ... , coming from 
n + 1 loop effects. Such an analysis was carried out in Ref. [15] for the leading log three-loop 
contribution. Since we have now clarified the short- distance two-loop behavior of aj^^, it is 
appropriate for us to revisit the issue of higher orders and refine the previous study. 

An interesting subtle feature that enters this RG analysis is the mixing of operators. We 
are interested in the OPE coefficient of the dimension 5 dipole operator jlaaplJ^ F"^. Leading 
logs contribute at the two-loop level due to QED corrections to the dipole operator (its 
anomalous dimension) as well as from two-loop mixing between the dipole operator and 
d — 6, four-fermion operators. A careful treatment of their mixing is important for the RG 
analysis. 

Before addressing the details of a|^^, it is useful to recall that a related QCD study was 
carried out about 25 years ago [27, 28] for the case of weak radiative decays which involve 
flavor changing gluomagnetic and electromagnetic dipole operators. Indeed, there QCD 
effects are very large and a RG analysis is essential. Later, because of the phenomenological 
importance of 6 — s'j, interest in such transition dipole operators increased, generating many 
studies and some controversies involving subtle issues regarding renormalization scheme 
dependence, 75 definition, operator set completeness etc. A brief discussion of those issues 
will provide guidance for our three-loop analysis. 




(^2) (P3) {P7) 




{F2) (F3) (F7) 

FIG. 5: Mixing of four-fermion operators V^/, A^f with H. The labeling of diagrams follows Ref. [31]. 

To discuss the issues that confronted radiative quark decays, we consider the two-loop 
Feynman diagrams in Fig. 5. These graphs describe mixing of the four-fermion operators, 

designated by Cg), with the dipole operator. Originally [28] only P2,3 and ^2,3 type diagrams 
were accounted for. Pj and Fj were later calculated in Refs. [29, 30] (those authors also 
accounted for one- loop mixing between gluomagnetic and electromagnetic dipole operators). 
Numerically, they did not cause much of an effect: about 7% in the mixing coefficients. The 
smallness of P7, F7 is related to the smallness of the internal fermion loop entering as a 
subgraph, which is basically a part of the photon vacuum polarization operator (it is also 
non- leading in a 1/Nc expansion). The smallness of this fermion loop was used both in 
[28] and [29, 30] to limit the number of four-fermion operators considered to a reduced set. 
The full set (originally introduced in [28]) contains penguin operators with right-handed 
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fermions that arise from left-handed ones due to the same fermion loop, see Fig. 6. So, their 
coefficients are also correspondingly small and could be neglected in early studies. 




(a) (b) 
FIG. 6: Renormalization of four-fermion operators: annihilation diagrams. 

The full operator basis was considered in later publications and we refer to [31] for a 
discussion of results and references to the literature. There, the renormalization scheme 
dependence and definition of 75 were shown to lead to different four-fermion operators. 
As explained in [31], to make the definition unambiguous, one has to redefine the four- 
fermion operators by adding to them dipole operators with appropriate coefficients. Those 
coefficients are fixed by the requirement that matrix elements of the redefined operators 
between fermion and fermion plus 7 states must vanish. In that way, consistency among 
different calculational approaches was restored.^ 

Here, wc note that the scheme independence can be understood in a simpler way. The 
basic point is that one-loop subgraphs for the two-loop diagrams in Fig. 5 arc finite and 
unambiguously fixed by the use of gauge Ward identities. A good example is the anomalous 
fermion triangle involving one axial-vector and two vector vertices where it is well known 
that the anomaly does not depend on the definition of 75. The logarithmic dependence on 
the normalization point (scale) which comes about due to the second loop integration is 
then clearly scheme independent. Below, we use this approach, originating from Ref. [28], 
to calculate all two-loop mixing. The results are consistent with those in Ref. [31] and with 
the explicit two-loop leading log calculations for given in section II. These consistencies 
provide a useful check on the analysis. 

Our calculation of the two and three loop leading logs differ, however, from the results 
in [15]. The disagreement can be traced to differences in the one and two loop anomalous 
dimension matrix elements. Details are given below. 

A. One- and two-loop results 

As we discussed in Sec. II the clectroweak contribution to the muon magnetic anomaly, 
fl/i = iOfi ~ 2)/2, can be represented as a sum over W, Z and Higgs bosons. In one-loop the 
Higgs contribution is negligible and aj^^(l-loop) given in Eq. (2) is a sum of a|j^)(l-loop) 



^ We disagree, however, with the statement in [31] about scheme dependence in case of the reduced set of 
four-fermion operators. In our view it is again related to the definition of four-fermion operators in the 
full set. 
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and a(f) (1-loop), 



af)(l-loop) = -^^ 
a(f)(l-loop) = -^^- 



10 



(68) 



where we denote the axial- vector and vector couplings of Z to the muon by g^^ and gy. In 
the Standard Model, 



5a 



2r = -1 



= 2ll - AQ,s 



w 



4s^ - 1 , 



(69) 



and gy is numerically very small. 

At the two-loop level electromagnetic corrections arc enhanced by log m^/m^. For a'^^\ 
the logarithmic part of the full two-loop result [14] is particularly simple, 



)(2-loop),, = -4 ^ In ^ ar(l-loop). 

TT Tfli, 



(70) 



As we discuss below, it just reflects the anomalous dimension of the corresponding dipole 
operator. In the case of a^^^ it is the only source of the logarithm at two loops. 

For a^^^ the situation is more complicated. In particular, Feynman diagrams without 
closed fermion loops give [10, 14] 



a^^-* (2- loop; no ferm. loops) 



87^2^/2 tt m 



-4 - In — a f Ul-loop) H ^^-^ • - In 



mz 



13 
9" 



IV) 



(71) 



In the second line we separated out the piece due to the anomalous dimension. 

We also have to add diagrams with closed fermion loops. The diagrams with the muon 
loops give 



ajf^ (2-loop; muon loops)^^ = 8^vt ' ^ ^ 



mz 



m,. 



-^N,{g>Xf - -N,{g^yf 



(72) 



where we introduced the factor equal to 1 for the muon loop just to distinguish between 
contributions with and without closed fermion loops. This generalizes calculations in [10, 
12, 13] by including the second term proportional to {gyY in Eq. (72). This term vanishes at 

arises from the induced coupling 



•'W 



1/4. In Eq. (72) the first term proportional to [gj^ 



of a Z with two photons via triangle diagrams, see the diagrams F2, F3 in Fig. 5. The 
second term, from the vector coupling, corresponds to the 7-Z mixing via a muon loop, see 
the diagram F-j in Fig. 5. 

Fermions other than muon contribute only via closed loops in two-loop order. Including 
their effect leads to a generalization of Eq. (72) to 



a|f^ (2-loop; ferm. loops)^^ 

^ Gf,ml a^^^ mz 
87r2-\/2 TT ^ {m/,m^} 



-Q9'A9iNfQ) + ^gWvN}Q} 



J 



(73) 
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where we introduced the notation 



{nif, m^} = max{m/, m^} . 



(74) 



Moreover, Qf is the electric charge of the fermion, Nf — 1 for leptons and Nf — Nc — 3 for 
quarks, and 



g{ = 27) , g(, = 27| - As'^Qf . 



(75) 



It is imphed that rrif ^ mz in Eq. (73); so, it does not include the top quark contribution 
which is part of the nonlogarithmic, NLL, terms. 

The closed loop contribution (73) and the last term in the second line of Eq. (71) are 
due to the two-loop mixing with four-fermion operators to be discussed below. Overall, the 
two- loop result for the sum of a^^^^ and a^^^ is (the form is slightly different than Eq. (7) in 
Sec. II, but equivalent) 



af'^)(2-loop) 



LL 



a 



In 



m,, 



+ E 

f=u,d,s,c,T,b 



Qg^ANfQ) + -9WvNfQf 



mz 



In ^ , (76) 



where we neglected the mass difference between W and Z (the ln(m^/mvy) terms are put 
into NLL contributions). The first term in Eq. (76) accounts for diagrams with muons 
and electrons and in the second term the sum is over all other fermions except top. The 
dependence on enters via gy^ gy . Our Eq. (76) differs somewhat from Eq. (25) in [15], 
as explained at the end of section IV C. 



B. Effective Lagrangian 

In the effective Lagrangian, represented as a sum over local operators normalized at a 
point fi, the anomalous magnetic moment of the fermion / is associated with the operator 
F'ajS f o""'^/ of dimension 5. Because of a chirality flip it enters with a coefficient proportional 
to the fermion mass m/, and it is convenient to include m/ in the definition of the operator, 

i^(/^) = -^[eF„/3/xa"'^/x]^, (77) 

where the electric charge e = y/ina is another factor included in the operator definition. 
Both the mass and the electric charge e arc /i dependent quantities in Eq. (77) but the 
running of the electric charge e is canceled by the wave function renormalization of the 
electromagnetic strength tensor Fap. The product eFa/B is RG invariant. 

The effective Lagrangian for flavor and parity preserving transitions can be written as 

J^M = Uf,)H{t,)+Y:cWi{A , (78) 
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where the second sum extends over d = Q four-fermion operators.^ The observable value of 
related to the coefficient h at the low normahzation point // = , 



,EW 



87r2V2 



(79) 



The relation (79) implies that only the operator H contributes to the matrix elements 
(/i|£efr(^/i)|At7) — a condition which we discussed above. 

The one-loop results (68) refer to the range of virtual momenta of order mw and mz, 
thus fixing the value of h at the high normalization point, 



10 



1 
3 



We choose the basis for the four-fermion operators defining them as 
The d = Q part of the effective Lagrangian can be written as 



E 



f,g,T=V,A 



(80) 



(81) 



(82) 



where Or-jg = Or-gf] so, the OPE coefficients c^'-^^ are symmetric under permutation of / 
and g . The tree-level Z exchange gives the initial data for the OPE coefficients. 



[r = v,A), 

where the vector and axial- vector couplings (^y^^ are given in Eq. (75). 

f f 



(83) 



FIG. 7: Contribution to from operators in C^r^{mw)- Both fermions / and /' must be charged 
so leptons do not contribute. 



In the case of W exchange the tree-level effective Lagrangian is 

V2 



uY{l - j5)ddj„{l - 75)m + {u^ c, d ^ s) 



(84) 



^ We omit here the d = 5 dipole operators for other fermions as well as chromomagnetic dipole operators 
for quarks, since they do not mix with H . 
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where we neglect CKM mixing. Such operators contribute to h only if all fermions are 
charged, so we can omit the leptonic part (see Fig. 7). One can use the Fierz transformation 
to put the operator in the following form: 



uY{l - 75)itrf7,,(l - 75)^ 



+2MtV(l - 75)McJr7,(l - 75)rf, 



(85) 



where are matrices of the color generators. The part with does not contribute to the 
magnetic moment (up to gluon corrections which we do not consider here) and neither does 
the parity breaking part, so C^^^ reduces to the form (82) with the initial data 



= C(m^) = ^ , (r = A) . 



(86) 



There are similar operators and coefficients for u, d substituted by c, s. These are the only 
operators that contribute to three loop leading log mixing in the W sector. 
The RG equations which allow us to calculate the running of /i(//) are 



d/x 2-n 



InHfJ') + E P^-J9 C^'^'(/^)^(/^ - rUfg) 



dc^;-^»(/i) _ _a(/i) 
d/x 27r 

da(/i) 



r',f',g' 



2tt 



(87) 
(88) 
(89) 



Here 7//, (^r-jg and Ir'^fg/ form the matrix of anomalous dimensions for operators H and 
Or-jg ■ That matrix is "block-triangular" : H does not mix with the d — 6 operators but 

the operators O^g do mix with H. The Pr-jg correspond to these mixings.^ The 9 functions 
in r.h.s. count only active fermions at the given jj,, with rufg denoting the maximal fermion 
mass in the operator Ofg , 

m/g = {m/, rrig} = max {m/, TTig} . (90) 

Perturbatively h — h^^^ + h^"^^ + h^"^^ + . . . where the raised index denotes the number 
of loops. In one-loop approximation h^^\fi) = h{M) where h(M) is given in Eq. (80). In 
two-loop order one can neglect in the r.h.s. of Eq. (87) the running of a, ch and c^''^^, and 
get 



a{M) 
2ti 



(91) 



Below we will compute and Pv-jg and then verify as a check that Eq. (91) matches the 
explicit two-loop calculations. We will then use the RG equations (87-89) to determine h 
in three loops. 



^ Note that our definition of anomalous dimensions differs from that in Refs. [15, 31] by a factor (— 1/2) . 
Also the normalization of four- fermion operators with f ^ g is different. 
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C. Anomalous dimensions and mixing of effective operators 



First, consider the anomalous dimension 7// of the dipole operator H . It is conveniently 
computed in the Landau gauge, with a photon propagator —i{g^u — k^k^/k'^)/k'^, since in 
this gauge there are no logs in Z factors. Thus, jh is given by the sum of three diagrams 
in Fig. 8 minus anomalous dimension 7^^ = 3 of the mass mju(/x) included into the definition 




(a) (b) (c) 

FIG. 8: Anomalous dimension of the dipole operator H. 



(77) of H. The result is 

7^ = -l-2-2-3 = -8, (92) 

where the numbers correspond to diagrams a, b, c and (— 7m)- 

More involved two-loop calculations are needed to determine the mixings Pr-jg of the 
four-operators (81) with H. The relevant diagrams are shown in Fig. 5. It is clear that the 
operator Or-jg mixes with H only when at least one of its fermionic indices coincides with 
the muonic one. 

Let us start with operator Oy-^^if with f ^ fi . It is the diagram F7 (plus, of course, a 
similar diagram where the virtual photon is coupled to the incoming muon leg) which defines 
• The fermion loop in this diagram (which is the same as in the photon polarization 
operator) produces 

where Q is the Euclidean momentum of the virtual photon. The factor in (93) cancels 
the photon propagator and for the second loop integration we get an expression similar to 
the one-loop Z boson exchange with the pure vector coupling, gy = I, gA = 0, up to the 
substitution 

"I -^-'^fQ>i^Ml. (94) 



m| + g2 247r2 
The ln(M^/(5^) can be represented as 

m 



= /" (95) 



with the range //^ ^ <^ M'^. That allows us to get the two-loop result from the one- loop 
one. From Eq. (80) we see that /i*^^) = 1/3 at gv — l,gA — 0. Thus, the diagram F7 produce 

e^NjQj 1 [M- dM' _ 2 e'NfQj ^ _ ^^^^ 



247r2 3V 9 I67r2 ^j? 
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in /i(^) . It gives for 

Pv;,f = g NfQf , / ^ , (97) 

where we accounted for another diagram similar to Fy . 

For the operator Oa;^/ with / 7^ // its mixing with H is given by the diagrams F2 and 
F3. The fermion loop in this case is the anomalous triangle with one axial- vector and two 
vector vertices. We need its kinematics when the momentum of the external photon tends 
to zero. The triangle then reduces to (see Eqs. (15) and (20)) 



e'NfQ 
27r2 



/ 



p _ Ml F _ Ml p 



(98) 



where q is the momentum of virtual photon and F^^^, = {\/2)e^i,^sF"^ is the dual of the 
external electromagnetic field F^ij. Integration over q in the second loop is logarithmic and 
produces 

(3A;,f = -^NjQ), f^^i. (99) 

For the pure muonic operators Oy-^h their mixing with H results from diagrams -^2,3,7 as 
well as diagrams i-*2,3,7 in Fig- 5. The F diagrams which are due to pairing of fermions from 
the same current in the four-fermion operators coincide (up to substitutions Nf — > A^^ = 1, 
Qf^Q^ = —l and the combinatorial factor two due to two ways of pairing) with those of 
the Or;/i/ operators, discussed above. 

We can transform the P diagrams, which are due to pairing of fermions from different 
currents, into F-type diagrams with closed fermionic loops using Fierz transformations, 

AYl5i^2 ^371^75^4 = ^ ^l7''^4 ^37i/^2 + ^ ^l7''75^4 ^37^75^2 

+'01'04 '03'02 - '0l75'04 '0375'02 ■ (100) 

We see that the P diagrams can be reduced to already calculated F ones (first two terms 
in the r.h.s. of Eqs. (100)) and to diagrams of the ^2^3 type where instead of products of 
axial-vector currents in the four-fermion operators we have scalar or pseudoscalar ones. 

Taken separately, the fermion triangles with the scalar and pseudoscalar vertices contain 
logs and produce double logs in the anomalous magnetic moment. But for the difference of 
scalar and pseudoscalar operators entering Eqs. (100) these double log terms cancel. What 
remains in this combination can be presented as a piece in the pseudoscalar triangle of the 
form 

-2^2<1^F^^- (101) 

The second loop integration is then simple. 

Altogether, it results in the following mixings of Ot;^ih with H 

4 / 4\ 22 8 

/3v,. = -g-6 + 4 + 27V,(--)=----iV,. 



4 94 
Pa-,, = -- - 6 - 4 + 27V^ (-6) - - y - 12A^^ ■ (102) 
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The numbers written after the first equahty signs display a decomposition in terms of closed 
loops: vector, axial-vector and scalar plus pseudoscalar loops. The latter piece has different 
signs for /5y;^^ and (^a-hh ■ We can unify the expressions (97), (99) and (102) as 

4 4 / 22\ 

Pa-j, = S>} ■ {-6N,Ql) + 5^ . {-6NfQ}) + 5^ 5^ ■ (-^) . (103) 

Now we are well prepared to compare the RG analysis with the explicit calculations of 
two- loop effects. For the W exchange the two-loop expression (91) reduces to the term 
with the anomalous dimension = —8 which matches the result (70). In the case of a^^\ 
inputting in Eq. (91) the initial data (80), (83) and the mixings (103) we observe at = 171^ 
full agreement with the sum of Eq. (71) and (73). 

The total two-loop result for a'j^^ + a^p given in Eq. (76) differs from that in Ref. [15] in 
the term {4:/9)gYgyNfQ f. In [15] it is multiplied by a factor {—Qf). This error originated 
in Eq. (23) of [15] for the two-loop mixing of the operators V"^/ with H. In accordance with 
the diagram Fy in Fig. 5, the factor Qj^ in that equation should be substituted with QfQ^ . 



D. The third loop effect 



To determine h at the three-loop level from Eq. (87) we have to account for the running 
of «(//) and c^'^^{ijl) up to the first loop. 



a{iJ,) = a{M) 



^ a{M) , M 

1 + 4^ E^/ In 



27r 



{H,mf} 



c^;-^^(/i) 



27r 



El ;/9 
^r';f'g' 



rj, r'-f 



V(M)ln^^ 



r',f',9' 



{^^,mf'g'} 



Using this expressions as well as the two-loop solution (91) for h{ii) we find. 



87r2 



— + ^bfLf 

"I'lj, f 

M 



+ E/5r;/,c^^^^(M) 
/,9,r 



ifg 



+ E/5r;/.7:^,c^'^'^'(M)Lj: 



-fg 
f'g' 



Here 



^ifg = 



In^ 
In^ 



M 



{mi,mfg,m^} 
M 



+ 29{mi - {rufg, m^}) In 



nil 



In 



M 



{mfg,m^} nil' 



ifg ^ ^ 

{mfg,mfg^,nif,} 



+ 2e(m/Y-{m/g,m^}) In 



iT^f'g' 



In 



M 



{mfg,m^} m/ 



(104) 
(105) 



(106) 



(107) 
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Additional input, needed for the calculation of c\j' , is the anomalous dimension matrix 
for four-fermion operators. The anomalous dimension matrix 7r /'g' is determined from the 
one-loop diagrams in Fig. 9 and Fig. 6. The calculations are relatively straightforward: the 



(a) (b) 
FIG. 9: Renormalization of four-fermion operators: without annihilation. 

diagram Fig. 9(b) produces no logs in the Landau gauge and can be omitted, the diagram 
Fig. 6(a) reduces to Fig. 6(b) by means of the Fierz transformations (100). The result for 
nonvanishing entries in "Yr-j'g' 



A;fg 
^V;f'g' 

^V;f'g' 



-6QfQgSj,5l,, 



;i08) 



Let us now use the expression (106) to calculate the third loop for the W exchange. In 
this case the effective Lagrangian (84) at = rriy/ contains only quark operators. They do 
not mix directly with H so the second term in (106) does not contribute. In the last term 
it is the mixing of quark operators with Oy^q^ which provides a nonvanishing contribution. 
Another contribution comes from the first term in (106), associated with the anomalous 
dimension of H. Using the initial data (86) for cF'^^^M) after some simple algebra we arrive 
at 



+ - 



a 



10 



8 In 



2 f^W 



10 40 aimw) , thw 

In 

3 3 



TT 



m, 



(109) 



m 



f 



32 mw 1 nT'W 

— In In 

81 mq rric 



32 , r, mw 
— In^ — 
81 mq 



where we introduced mq as an effective IR cutoff for the light quark loops, q = u,d, s, 
implying mq is larger then We also added the first and second loop to make easier to 

follow a numerical comparison. Taking mq — 0.3 GeV, mc — 1.5 GeV, = 4.5 GeV we get 



^ Note that in sections II and III we used nis = 0.5 GeV different from m„ = nid = 0.3 GeV but here we 
make the simphfication ruu = ma = mg = mq. That difference has no numerical impact. 
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numerically, 



(W) 



a(mw) , 
1 - 26.5 ^ ^' + — (352 + 359 - 6) 



= [ 1 - 0.067 + 0.0038] , (110) 

showing that the third loop effect is quite small. It is dominated by the the anomalous 
dimension term oc '-fj^ (the first number 352 in term) and by the cross term oc ^nb 
between the anomalous dimension and running of a (the number 359), the four-fermion 
operators contribute very little, (—6). Moreover, the effect of the third loop becomes even 
smaller if we use «(m^) instead of a{mw) in the second loop: this changes the ^ynb term in 
the third loop, 359 -252 . 

Finally, we note that if we shift to the usual fine structure constant, a = 1/137.036, via 
the full leading log relation 

aimw) = « + ^ ^ NfQ} In ^ ~ 1/129, (111) 
37r ^ •' ruf 

it generates from the 0{a{mw)) terms in Eqs. (109) and (110) additional 0{a'^) contri- 
butions that amazingly cancel (numerically) the 0.0038 in Eq. (110). So, in terms of the 
expansion parameter a, the leading log three-loop corrections are essentially zero. Such a 
complete cancellation appears to be a numerical coincidence. 

The algebra is more tedious in the case of the Z exchange. We set = 1/4 which is a 
very good approximation numerically. It simplifies the analytic result due to vanishing of 
the leptonic vector couplings gy'^'^ at this value of s^, leaving us with fewer operators. We 
also combine the W and Z exchanges neglecting the mw, fnz mass difference — again quite 
a good approximation. (The difference will be in the three-loop NLL.) 

We present the final result for the three-loop part of in the form similar to that used 
in Ref. [15], 

aJ^'^(3-loop),, = ^ . ^ ^ (A, + A, + + B,) , (112) 
where Ai^g come from lepton and quark terms in Eq. (106) containing jjj, (3j 

= ^ In^ !^ - ^ In^ + In ^ In , (113) 
90 45 rriT- 5 m,- 



2662, oTTiz 11216, 2 1964, o "^z 24, mz , mz 

m 1 m 1 m 1 m m 

1215 mi, 1215 rric 405 uiq 5 



96 , mz , mz 48 , mz , mz , 32 mz , mz , 32 mz , mz 

m m m m 1 m m 1 m m , 

5 mc m^ 5 mq m^ 405 mc 135 mq 

and Bi 2 are due to the 7ij6, (3b terms involving running of a , 
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9 m„ 



mq 



mr 



mr 
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The B2 term can be fully absorbed into the two-loop part of given in Eq. (76) if a there 
is substituted by a{m^) instead of a{mz) which we used in our derivation above. 

Comparing with the results in Ref. [15] we see that our Bi coincides with their B but the 
sum Ai + Aq is somewhat different from A in [15]. This is due to a few reasons. One was 
already discussed: in Eq. (23) of [15] for the two-loop mixing of the operators T/^/ with H 
the factor Q'j should be changed to QfQ^ ■ In addition, there is a difference in the one-loop 
anomalous dimensions of four-fermion operators. In [15] an extra factor 2 is ascribed to the 
penguin diagrams in Fig. 6. To correct this the factor 1/2 should be introduced in Eqs. (33), 
(36), (38-41) of [15] and in Eq. (35) the 52/3 should be changed to 44/3. 

Numerically, 



1696, 



4 =-507, 5i = -774, 



B2 = 1916. 



(115) 



where we used the same values for the quark masses as above. Altogether, the three-loop 
correction is 

a^W(3-loop) = a^W(l-loop) (^^^ (A + Aq + B,) ~ 0.4 x 10"^^ , (116) 

where it is implied that a{m^) is used for the two-loop part. The numerical value is close to 
that given in [15]. The reason for this final agreement is that the A-part of the three-loop 
result is numerically dominated by 7!^ and large mixing of axial operators with the dipole 
operator, followed by the running of the dipole. These pieces, as well as the S-part, are 
correct in [15]. If a{mz) is used in the two-loop part the third loop is somewhat larger, 
af^{3-loop) ~ 2.4 X 10"". 

Of course, in both cases one must reevaluate a|^^(2-loop) with the shifted coupling at 
scale m^ or mz- In that way, scale insensitivity is restored. In addition, because the effective 
couplings are larger than the usual fine structure constant, a, a transition to this a in the 
two-loop part of induces additional negative contributions. Remarkably, those negative 
contributions cancel with the above explicit three-loop results to about 0.1 x 10~^^. (The 
cancellation is similar and of course related to the even more complete cancellation pointed 
out for the W contribution alone.) Hence, to a good approximation, the leading-log higher 
order contribution is zero or at least negligible. 
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V. SUMMARY 



Having addressed a variety of computational issues, 



Small, previously neglected, 2-loop contributions suppressed by factors of (1 — 4s^) 
that come from 7-Z mixing and the renormalization of 6w, 



• Strong interaction modifications of quark loop diagrams, and 

• Leading log 3-loop effects, 

we are now in a position to update the Standard Model prediction for a^"^ and assess its 
degree of uncertainty. 

Small effects due to 7-Z mixing and our choice of 9w renormalization have now been 
included in Eqs. (7) and (8). Because of the (1 — 4s^) suppression factor, their total impact 
is rather small, shifting the value of down by about 0.4 x 10^^^. 

More important are strong interaction effects on the quark triangle diagrams in Fig. 3, 
particularly in the case of light quarks. It was shown that short distance contributions are 
unmodified (thereby, hopefully, eliminating controversy in the literature). However, QCD 
can affect their long-distance properties. In the case of the first generation of fermions a 
detailed operator product expansion analysis and effective field theory calculation led to a 
shift relative to the free quark calculation (with constituent quark mass) by 

AaJW[e, u, d]QCB - ^af^[e, u, djfree quacks = +1.7 X 10-'\ (117) 

For the second generation, comparison of the free quark calculation with the more precise 
evaluation in Eq. (65) shows no significant numerical difference. However, the more refined 
analysis now indicates very little theoretical uncertainty. So, the total hadronic uncertainties 
in would seem to be well covered by an uncertainty of ±1 x 10^^^ or even less. 

Finally, after a detailed renormalization group analysis, the leading log three-loop contri- 
butions turned out to be extremely small. In fact, they are consistent with zero, to our level 
of accuracy ~ 0.1 x 10~^^, due to a remarkable cancellation between anomalous dimensions 
and running coupling effects. Uncalculated three-loop NLL contributions are expected to 
be of order 

''""J^^^yto^^g.lO-", (118) 



which is negligible unless enhanced by an enormous factor. Nevertheless, we assign an overall 
uncertainty of ±0.2 x 10^^^ to a^"^ for uncalculated three-loop NLL contributions. 

So, in total we find a small shift in a^"^ (for itlh — 150 GeV) from the previously quoted 
value of 152(4) x 10~^^ to a slightly larger (but consistent) value 

= 154(1)(2) X 10-^^ (119) 

where the first error corresponds to hadronic loop uncertainties and the second to an allowed 
Higgs mass range of 114 GeV ^ mn ^ 250 GeV, the current top mass uncertainty and 
unknown three-loop effects. 
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